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Abstract 

Spatial correlations in sheared isothermal liquids for both elastic and granular cases are theoreti- 
cally investigated. Using the generalized fluctuating hydrodynamics, correlation functions for both 
the microscopic scale and the macroscopic scale are obtained. We find the existence of long-range 
correlations obeying power laws. The validity of our theoretical predictions has been verifled from 
the molecular dynamics simulation. 
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I. INTRODUCTION 



Liquids consist of assemblies of many particles. There is long history to study molecular 



liquids which consist of non-dissipative simple molecules 12|. In these days, there has 



been rapid growing interest in granular liquids which are made of granular assemb 



dissipative macroscopic particles 3, ^, Q, y, y, [lo|, [ll|, [l2, 14, Q, 



16 



les, I.e., 



17 



18|. 



Although a molecular liquid is fluctuated around an equilibrium state, the absence of the 
equilibrium state of granular assemblies makes granular materials unlike usual materials [l^ . 
Nevertheless, we have recognized that statistical properties and hydrodynamic behaviors of 



moderate dense and nearly elastic 



analysis of the kinetic theory 



granu 



20 



21 



ar g ases 



22 



23 



are considerably understood from the 



24|. 



Boltzmann-Enskog theory is often used in describing moderate dense gases including 



granular assemblies, where the assumption of molecular chaos is used 



21 



22 



24S\ . In spite of the success of qualitative description of hydrodynamic behaviors based on 



Boltzmann-Enskog theory, significant roles of correlations have been recognized in these 
days. Indeed, it is well known that correlated collisions cause significant differences in 
constitutive equations even in equilibrium gases. In particular, recently a number of papers 



on long-time tails in current correlation functions for granu 



which are directly related to the transport coefficients 



11 



ar liquids have been published, 



12 



15 



3,0, 



25 



26 



27j. On 



ong-ra nge correlations in sheared elastic fluids 



33 



34 



35 



36, 



371] . The existence of a similar 



;he other hand, we know the existence of 

28, 2^, 30 1 and heat conduction systems 31 
long-range correlation has been observed even in a simulation of sheared granular fluids 38|. 
In spite of the indication of the existence of the long-range correlations in sheared fluids, 
we do not have any consistent theory in describing a structure factor for both the particle 
scale and the hydrodynamic scale. Indeed, a theoretical prediction of spatial correlations in 
the molecular scale for relatively dense liquids under a shear is not consistent with the long- 
range correlations [sQ^. Thus, we still do not understand the details of spatial correlations 
even in the case of molecular liquids. 



We often use the mode-coupling theory (MOT) in describing dense hquids [40|, |4l|, |42 |. 
Even when we are interested in sheared dense granular liquids, the framework of MOT can be 
used [17]. However, MOT is not a closed theory, because it needs to determine the structure 
factor by another method. Thus, to obtain the structure factor or the pair-correlation 
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function is an important issue for the description of dense granular liquids. 

The purpose of this paper is, thus, to clarify the spatial correlation functions for sheared 
isothermal liquids. Based on the generalized fluctuating hydrodynamics 43|, |44|, we will 
demonstrate that the long-range correlation obeying a power law is consistent with the 
short-range structure obtained by a liquid theory. In the next section, we will summarize the 
outline of the generalized fluctuating hydrodynamics. In Section III, we will show properties 
of a set of linearized equations around a uniform shear flow (USF) based on the generalized 
fluctuating hydrodynamics. In Section IV, we will calculate the spatial correlation functions 
and their asymptotic forms in the long-range limit. In Section V, we will compare our 
results with those of the molecular dynamics simulation. In Section VI, we will discuss 
and conclude our results. In Appendix A, we summarize the form of the pair-correlation 
function go{r,e) for homogeneous unsheared states obtained by Lutsko [gI]. In Appendix B, 
we briefly summarize the transformation between the Cartesian coordinate and the oblique 
coordinate. In Appendix C, we present the explicit form of which is the Fourier 

component of the density correlation. In Appendix D, we show the explicit form of Cpp{k) 
which is the Fourier component of the momentum correlation. In Appendix E, we evaluate 
the asymptotic form of a function Aj(r) which is needed for the calculation of correlation 
functions. 

II. GENERALIZED FLUCTUATING HYDRODYNAMICS 

We consider three-dimensional systems consisting of identical smooth and hard spher- 
ical particles confined in the volume V = under a shear flow with the shear rate 7, where 
each particle has the mass m and the diameter a. The particles collide instantaneously with 
each other by a restitution constant e which is less than unity for granular particles, and 
is equal to unity for molecular liquids. Let us assume that the restitution coefficient e is a 
constant. 



The spatial corre 



dynamics 



28|, 



29 



30 



ations in sheared fluids are investigated by using the fluctuating hydro- 



38l | which can be used for the description of hydrodynamic behaviors. 



However, it is possible to generalize the fluctuating hydrodynamics which can cover the scale 



around the particle diameter. The generalized hydrodynamic equations 43|, |4J, |45|, |46|] for 
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isothermal liquids characterized by a uniform temperature T are given by 



dtn + V ■ (nu) = 0, (1) 
dtUa + u^VpUp + -Va/x + —V p{T.^p + Ef^) = 0, (2) 

where n and Ua are the number density and a-component of the velocity field, respectively. 
Here, we have introduced the generalized chemical potential or the effective pressure 



lnn-ldr'Cir-r',e,^)5nir',t) + 



(3) 



where C{r — r',e, 7) is the two-particle direct correlation function, which satisfies 
nC(A;,e,7) = 1 — S{k,e,'y)~^ with the structure factor S{k,e,'y). E^^(r,i) is the viscous 
stress tensor given by 

^apir,t) ^ - J dr'[r){r-r',e)i^eap{r',t) + epa{r',t)-'^e^j{r',t)Sai3^ 

-C{r-r',e)e^^{r',t)Sap], (4) 

where iapir^t) = {'S/aUf3{r,t) + V pUair^t)} /2. T,^p{r,t) is the random part of the stress 
tensor satisfying (E^^(r,t)) = 0, and 

(Ef^(r, m^sir', t')) = 2TS{t - t'){rj{r - r', e)A^p^s + C{r - r', e)S^f,S,s} (5) 

with Aaissy = SajSf^s + ^aS^pj — ^Sa/sS^s/^- Here, we have used Einstein's sum rule on the 
Greek subscript. The generahzed shear viscosity rj{r, e) and the generahzed bulk viscosity 
({r, e) are represented by i'*{k, e), f|(A;, e) and Enskog's mean free time 

with the radial distribution function go{a, e) at contact as i^i (/c, e) = {m,noaH^^)~^{({k, e) + 
4:r]{k, e)/3), and (/c, e) = {mnQ(TH'^)~^r]{k, e), where no, r]{k, e) and C,{k, e) are respectively 
the average number density, Fourier transforms of r){r,e) and ({r,e). It is known that 

u^kje) and i'2{k,e) are respectively given by 

i^l{k,l) = 2{l-Mk) + 2Mk))/{3k'), (7) 

and 

u*{k,l) = 2{l-jo{k)-h{k))/{3e) (8) 



for elastic hard spherical particles, where ji{k) with / = or 2 is the l-th. order spherical 



Bessel function ka, 



44 



46 



47 



481]. Although we do not know how z/J(A;,e) and z/2(fc,e) 



depend on e, the explicit e-dependences of h'l{k,e) and h'2{k,e) are not important in this 
paper. Therefore, we are keeping discussion without their explicit forms. 

It should be noted that sheared corrections to the structure factor or the pair-correlation 
function can be obtained within this theoretical framework, though the equilibrium or the un- 
sheared structure factor So{k, e) = S{k, e, 7 = 0) should be determined by another method. 
We adopt an approximate expression of the pair-correlation function for unsheared granular 
liquids obtained by Lutsko which covers the equilibrium pair-correlation in the elastic 
limit (see Appendix A). 

This set of equations ([I])-® is a reasonable starting point for isothermal molecular liq- 
uids (e = 1), once we use appropriate generalized transport coefficients and So{k,e = 1). 
However, the validity of this set of equations in describing isothermal sheared granular fluids 
might be controversial. Indeed, nobody has used the generalized fluctuating hydrodynamics 
for granular liquids, because we know that fluctuations cannot be characterized by Gaussian. 
In addition, there are no explicit characteristics of granular liquids in the set of equations 
except for e-dependence of 5'(/c,e,7) and the generalized transport coefficients. Even when 
we use the hydrodynamic equations for granular liquids, we should consider an equation for 
the granular temperature. 

Let us answer the above critical points to validate the generalized fluctuating hydro- 
dynamics in eqs. ([I])-©. First, the granular temperature T is not a true hydrodynamic 
variable but a relatively fast variable because of the collisional energy loss in granular sys- 
tems. Thus, we expect the fast relaxation to a steady state of the temperature for sheared 
granular liquids, which satisfies T oc 7^/(1 — e^). Second, an isothermal situation is easily 
realized by the balance between the viscous heating and the collisional energy loss. In partic- 
ular, it is known that USF is stable for small and nearly elastic systems under Lees-Edwards 
boundary condition. Even in physical situations, the heat conduction is not important in 
the bulk region far from the boundary. In these situations, we may assume that sheared 
granular liquids are nearly isothermal. We also indicate that our previous studies clarify 
the formal similarities between sheared granular liquids and sheared molecular liquids at a 



constant temperature 



16 



171] . Through previous studies, we have recognized that the most 
important issue is to determine S{k, e, 7) for granular liquids [itI, which can be determined 
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within the framework of the generahzed fluctuating hydrodynamics. Third, the fluctuating 
lydrodynamics has been used in describing granular hydrodynamics for freely cooling cases 



49|. Thus, we believe that the set of equations ([I])-® can be used even for sheared gran- 



ular liquids. Although there are uncovered regions of our approach in the description of 
granular liquids, we expect that our approach can capture some aspects of sheared granular 
liquids. The validity of the model will be tested from the comparison between our theoretical 
prediction and the direct simulation of granular assemblies. 



III. LINEARIZED EQUATIONS AROUND UNIFORM SHEAR FLOW AND 
THEIR SOLUTIONS 

In this section, we analyze a set of the linearized equations around USF. As mentioned 
in the previous section, we assume that USF is stable. Thus, we only need to solve the 
linearized equations. This section consists of two parts. In the first part, we summarize the 
expression of the linearized equations. In the second part, we explicitly write the solutions 
of the linearized equations. 



A. Linearized equations 

Let us introduce the fluctuations of the hydrodynamic fields 6n{r,t) = n{r,t) — no, 
6u{r,t) = u{r,t) — c{r,t) with Ca{r) = 'jySa^x, and the non-dimensionalized vector z{r,t), 
whose Fourier transform is given by 

z^{q,t) = {6n{q,t),6u,{q,t)/{t-^'a'),6uy{q,t)/{t-E'a''),6u,{q,t)/{t~^'a')). (9) 

Here, the Greek suffix a denotes the Cartesian component. From eqs. ([I]), ([2]), and we 
obtain the linearized evolution equation for z{k,t) = z{q,t) as 

{dt-rk.^Yz + l-~z = R, (10) 

where the time, the wavenumber and the shear rate have been non-dimensionalized by 
t = tEt, q = k/cr, and 7 = Y I'^e, respectively. Here, the matrix L is expanded as 

L = Lo + 7*Li + ---, (11) 
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where Lq and Li are respectively given by 



Lo 



p*iky 



[yl ^2)kzkx 



rioa^iky 



(z/* i'2^kxky 



1^2 ) kz ky 



UQa^ikz 
{ul - vl)kxkz 
{ul - ul)kykz 
[ul - ul)kl + ulk^ 



[12) 



Li 
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(13) 



with p\ = p*{k,e) = ASQ{k,e) ^ and A = T j {mriQa^t^). The random vector R has four 
components 



Ri = 0, 

Rc+i = {mnoaH^'^)'''^iakpE^Jk,i), 



(14) 



where a = 1, 2, 3 respectively correspond to x, y and z. Although 7 dependence of S{k, e, 7) 
should appear in Li, we simply ignore such terms. The validity of this simplification will be 
checked from the comparison of the results with our simulation. 



B. Solution of the linearized equations 



its calculation can be 



It is straightforward to solve eq. fllOl) . As mentioned in ref. 
simplified if we introduce the transformation from the Cartesian coordinate to the oblique 
coordinate to decompose the longitudinal modes and the transverse mode. However, the 
obtained results in the oblique coordinate are rather confusing, because of too many suffices, 
n addition, the results are basically the same as those obtained by Lutsko and Dufty 



28l | with replacing the transport coefficients by the generalized transport coefficients, and 



ignoring terms related with the fluctuation of the temperature. In this paper, thus, we 
only present the final results in the Cartesian coordinate. The transformation between the 
Cartesian coordinate and the oblique coordinate is explained in Appendix B. 
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The solution of eq. (fTOl) can be formally represented by 

4 rt 



7=1 -^-"^ 



(15) 



where 



i^^^\k,t} = ^p^^\k)exp[- [ rff A(^)(fc(7*f ))], (16) 

Jo 

F'^^\k,t) = ip'^^\k) ■ R{k,t) (17) 

with k{f) = {kx,ky + fkrc,kz). Here, we have introduced the right eigenvectors %j)^'-'\k), 
the associated biorthogonal vectors, i.e. the left eigenvectors (p^^\k), and the eigenvalues 
\'^^\k) satisfying 

(^-l^*k^-^ + • V^^'Hfc) = \^^\k)il^'^^\k). (18) 

We also note that ^^^^(fc) and c^(-')(fc) satisfy ip^'^k) ■ ip^^^k) = 6ij. 

In order to obtain ■^'•-'^(fc), (fi^^\k), and X^^\k), we use the expansions 

V^o-) = 4^') + Yi^? + ■■■ , (19) 

= (^(^■)+7V? + ---, (20) 
A^-) = Ai^VrAl'V-- - , (21) 

in terms of 7*. 

We should note that the perturbation in terms of 7* is not the expansion from an un- 
sheared state of granular liquids. Indeed, it is well-known that properties of sheared granular 
liquids completely differ from those of freely cooling granular liquids. In the case of sheared 
granular liquids, we obtain the relation 7* ~ 'j/VT ~ Vl — from the balance between 
the viscous heating and the collisional energy loss. Thus, the expansion in terms of 7* can 



be regarded as that by small inelasticity 501 ]. 



Substituting eqs. (fT9|) . (l20l) and (12T|) into eq. (fT8|) . we obtain the zeroth and the first 
order perturbations as 



;Lo-A^^^l)-<^(fc) = 0, (22) 

- -^ I 



(L(o) - A^^-^) ■ + f -lk.4r + Li - ) . il^ll\k) = 0. (23) 
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Solving these equations, we obtain the eigenvalues 





\ _l_ A,* ^ 
= ^+ + 7 ^2 




A{2) 




'<^'\k), 


A{3) 


= e)A;^ — 


k k 


A(4) 







(24) 

(25) 

(26) 
(27) 



within the approximation up to 0(7*), where we have introduced 



iy*{k, e)k' + ^{ul{k,e)k^)^ - Anoa^p*{k,e)k^ 
= 2 ' ^ ^ 

_ iy*{k,e)P - ^{ul{k,e)k^y - 4:noa^p*{k,e)k'' 
^- — 2 ' ^ ^ 

e\k) ^ ^ + ^kd,p*{k, e), C'^\k) ^ ^ + ^kd,p*ik, e), (30) 



and 



iV+ = -noaY{k, e)k'^ + A^, Nl = -noa^p*{k, e)k'^ + (31) 
Similarly, we obtain the right eigenvectors 

^<^)^ = i-(,,„„.3,A4,A4.A4). (33) 

^(3) = *(3)+M(fc)*(^), (34) 

^/;(^) = (35) 



(36) 



and the left eigenvectors 





= ]^ (^iA;p*(A;,e),A+^,A4 


ky . 


kA 

k) 




= -^[ikp*{k,e),\J^,\_ 




kA 
k) 


<^(3) 










= -M(fc)#(^) + *(^), 







(37) 

(38) 
(39) 
(40) 
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where we have used 



and 



with k± = k"^ — ky. 



VI>(3)T ^ ^(3) ^ fr. ^ _kyK\ f.^. 

\' kk^' k' kk± J ' ^ ^ 

(o.|i.O.-^). (42) 



kk 

M{k) = -^^tan-\ky/k±) (43) 



IV. CORRELATION FUNCTIONS 

This section is the main part of this paper, in which we present the exphcit forms of 
correlation functions. This section consists of two parts. The first part summarizes the 
general results of correlation functions. In the second part, we evaluate integrals included 
in correlation functions to extract the long-range behaviors of correlations. 

A. General results for correlation functions 

Let us introduce the correlation functions Cii{k,t) which satisfy 

{S,{k,t)zi{k',t)) = i27if6\k + k')Q,{k,i). (44) 

Note that we do not use Einstein's sum rule for Latin subscripts. Substituting eq. f|T5l) into 
this equation, we obtain 

POO _ 

Cu{k,t)= / dtJ2i'f\k,t)i>t\-k,t)F^'"'\kirt)), (45) 

where ipf\k,i) is the i-th. component of -j/?*" \k,i), and F^''''^\k) satisfies 

(F«(fc,r)F(™)(fc',t')> = {27rfS''{k + k')6{t - F)F^'^\k). (46) 
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From eqs. flTll) . f[T7|) and fl57l) - fllDl) . the explicit forms of are given by 



= F(^i)(fc,t) = -2AfcV*(A:,e)^, 
= 2AA;V2*(A;,e), 

F^^^\k) = =M(fe)F(=^3)(fc). (47) 

Thus, we can calculate any spatial correlation functions. 
Let us explicitly write the spatial density correlation 

Cnn{r,t) = {Sn{r + r',i)Sn{r',t))=a-'^Cn{r,t) (48) 

and the spatial momentum correlation 

Cpp(r, i) = {p{r + r', i) ■ p{r\ t)) ~ {mn^a ltEf{C22{k, t) + Cgslfc, t) + C^^{k, i)] (49) 

with the momentum density p{r,i) = mn{r,t)5u{r,t). In eq. (l49l) . we adopt the approxi- 
mation p{r,i) ~ mnHSu{r,t). With the help of the inverse Fourier transform 

Cnn{r,t) = J ^C'„„(fe,^)e-^''■^ (50) 

Cpp{r,i) = J j^^^Cpp{k,i)e (51) 

and eq. fH5l) . we obtain the steady solutions (see Appendices C and D), 

Cnnik) = no{So(A;,e)+7*Ai(fc)}, (52) 
Cpp{k) = [mnoa/tEYi-D.ik) + D,{k) - D,{k)}, (53) 

where C„,„(A;) = limt^ao Cnn{k,i), and So{k,e) is the structure factor for unsheared case. 
Dj{k) with j = 2, 3, 4 in eq. ( l53l) are given by 

b^ik) = -A|l + 7*A2(fc)}, (54) 
D,{k) = A|l + 27*A3(fc)}, (55) 
L)4(fc) = + 27*A4(fc)} . (56) 
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The derivation of eqs. (I52l) - (l56l) and the exphcit forms of Aj{k,t) with j = 1,2,3,4 are 
presented in Appendix O and [Dl 
From (!50l) -( l56l) . we finally obtain 



Cnnir) = n^{(?o(r,e)-l + 7*Ai(r)} + no5(r), (57) 
Cppir) = ^ [35(r) + a~^* {/^^{r) + /\^{r) + /\,{r)}] , (58) 

TflTlQ 

where go{r,e) is the pair-correlation function for unsheared cases, and 

A.(r) = I ^M^y''"'"- (59) 

The determination of gQ{r,e) or SQ{k,e) for e < 1 is highly nontrivial. Indeed, we can- 
not keep any homogeneous cooling state (HCS) without artificial controls of the systems. 
However, Lutsko 6| obtained an approximate expression of go{r, e) for HCS. He also verified 
that his approximate expression works well from the comparison between the theory and the 
simulation of HCS. As stated in Section II, we adopt his expression in this paper (Appendix 
A). We also note that the structure of the liquids in eq. (l57j) can be represented by the 
linear contribution of the homogeneous terms and the sheared term. 

B. Long-range correlation 

Let us demonstrate the existence of the long-range correlation in C„„,(r) and Cpp{r). Let 
the angular average of any function /(r) be denoted by /(r) = J dQf{r)/{ATT). As shown 
in Appendix E, the asymptotic forms of Aj(r) (j = 1, 2, 3) satisfy 

Ai(r) oc r-"/^ r > 1^, (60) 

A2(r) oc r-"/^ r > (61) 

Asir) oc r~^/^ r > 1^, (62) 

A4(r) oc r^^/^ r > /c, (63) 

where Ic = a j ^pf" . Substituting these results into (1571) and ( |58l) the long-range parts of 
Cnnij) and C'pp(r) respectively satisfy 

-11/3 



Cnn{r) « ( f ) , > /c, (64) 



Cvv{"^) [j\ ' r-^U. (65) 
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Cnn(r) 



n 



0.6 




FIG. 1: The density correlation function Cnn{f) for the volume fraction (p = 0.185 with L = 89(7 
and e = 0.83 as a function of the distance r. The solid line represents go{r, e) obtained by Lutsko[6t]. 



This long-range correlation Cpp{r) is known for isotherma 
and has been verified in sheared dilute granular fiuids [38 1 



sheared elastic fiuids |28|, l29|, l30| , 



V. COMPARISON BETWEEN THEORY AND SIMULATION 

To verify the validity of our theoretical prediction, we perform the event-driven molecular 
dynamics simulation for three-dimensional hard spheres. In our simulation the time scale 
is mesured by tq = a^Jm/TQ where Tq is the averaged initial temperature. Particles are 
confined in a cell under the Lees-Edwards boundary condition, in which each linear dimension 
is L. The number of particles is not fixed in our simulation, but we control the system size 
L and the volume fraction cj) as well as the restitution constant e. The initial state at the 
time t = is the equilibrium state, and we will show the correlation functions at t = 20ro 
for = 0.50 and 0.37, and t = 40ro for = 0.185 as Cnnir) and Cpp{r), where the system is 
considered in a steady state. We also choose the shear rate to keep the steady temperature 
unity in the dimensionless unit, except for the data of Fig. HI where the temperature in the 
steady state is I.OTq for 7 = 0.92tq^ or A.OTq for 7 = 1.84ro"^ 

Figures [IH3] show the behaviors of C„„(r) for r < 5cr at = 0.185,0.37, and 0.50, re- 
spectively. The restitution coefficient e is 0.83 for Fig. [H and e is 0.90 for Figs. [2] and [31 
The solid line represents go{r, e) obtained by Lutsko 6| without any fitting parameters. Al- 
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Cnn(r) 




1 1.5 2 2.5 3 3.5 4 4.5 5 



r/o 



FIG. 2: The density correlation function C„„(r) for the volume fraction (j) = 0.37 with L = 72a 
and e = 0.90 as a function of the distance r. The solid line expresses go{r, e) obtained by Lutsko[6t]. 



Cn„(r) 



n, 




r/CT 



FIG. 3: The density correlation function Cnnif) for the volume fraction (p = 0.50 with L = 32.5(7 
and e = 0.90 as a function of the distance r. The solid line expresses go(^) e) obtained by Lutsko 6|]. 

though we omit contributions of shear rate to ^^^(r) because of its simplicity, the agreement 
between the results of our simulation and our theory seems to be perfect. Thus, it is hard 
to find any contributions of the shear in the short-range structure of the density correlation 
function, as is known in dense elastic liquids. 

However, the above results do not mean that contributions of the shear to the density 
correlation function are not important. Indeed, we find the existence of a power law tail 



14 



Cnn(r) 



0.01 



0.0001 



■11/3 



■ QoM)-! 
7=0.92 To 




\ 



\^ = y = 1 .84 To 

^\ 
„ ■ \ 



10 

r/a 



100 



FIG. 4: The double-log plot of Cnnif) for the volume fraction = 0.185 with L = 89o", and 
e = 0.83 as a function of the distance r. Here, the solid line is that for homogeneous case obtained 
by Lutskojo^, and the plotted data are obtained in the case of 7 = 0.92t-^ and 1.84rQ-^ 

which might be consistent with the theoretical prediction C'„„(r) ~ /--n/s [y^ Qq f[64l) for 
r > (T (see Fig. H for = 0.185, and 7 = 0.92ro"^ and IMt^^). It should be noted 
that the shear rate dependence of Cnn{r) cannot be observed, because the scaled shear rate 
7* ~ i/Vr ~ Vl — in eq. (l57j) is independent of 7 in the steady state. The range of 
the tail obeying a power law is not wide enough to verify the theoretical prediction in Fig. 
m because of the large statistical errors. However, as will be shown in Figs. [7] and 8, the 
numerical data for Cpp{r) are consistent with the theoretical prediction. To confirm the 
quantitative accuracy of our theory in (7„„(r), we need more extensive simulations to reduce 
the statistical errors. 

Figures [5] and E] shows the behavior of Cnnij) for = 0.37 and 0.50, respectively. We use 
e = 0.90 and L = 72a for = 0.37, and e = 0.90 and L = 32.5a for = 0.50. Although we 
can identify the existence of a power law tail, it is hard to verify whether the tail satisfies 
fQj^ — Q_37 Q_5g because of the large statistical errors. 

The long-range momentum correlation function which satisfies r~^/^ can be observed 
more clearly than the case of the density correlation function as shown in Figs. [7| and [HI 
Figure His the result of our simulation for (p = 0.37 and e = 0.90. Figure [H] shows the result 
for = 0.50 and e = 0.90. Although there is an apparent finite size effect, Cpp{r) can be 
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FIG. 5: The double-log plot of Cnn{f) for the volume fraction (j) = 0.37 with L = 72a, e = 0.90 as a 
function of the distance r. Here, the solid line is that for homogeneous case obtained by Lutsko. 
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FIG. 6: The double-log plot of C„„(r) for the volume fraction (j) = 0.50 with L = 32. 5o", e = 0.90 
as a function of the distance r. Here, the solid line is that for homogeneous case obtained by 



Lutsko. 



an universal function of r/L, and decays faster than power-law function for r > 0.4L (see 
Figs. 7 and 8). We will have to check whether the oscillation of Cpp{r) in small r can be 
understood by our theory. To obtain the complete forms of Cnn{r) and C*pp(r) we need to 
solve the eigenvalue problem of inelastic Enskog operator and determine I'Kk, e) and ^2(^5 ^) 



16 



(mToG-'') 



0.1 



0.01 



0.001 



0.0001 



1e-05 



■5/3 



1 




100 



FIG. 7: The double-log plot of Cpp{r) for the volume fraction ^ = 0.37 with L/a = 36,72 and 
e = 0.90 as a function of the distance r. Here, a is a fitting parameter, which is 0.4 for L = 36(7 or 
1.0 for L = 72a. 
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FIG. 8: The double-log plot of Cpp{r) for the volume fraction cf) = 0.50 with L = 32.5(7 and e = 0.90 
as a function of the distance r. 

explicitly. This will be our future work. 

In this paper, we mainly focus on the results of sheared granular liquids, because gran- 
ular liquids are more non-trivial than sheared isothermal elastic liquids, in which particles 
collide with each other without any loss of energy. We also perform the molecular dynamics 
simulation for sheared elastic liquids with the velocity scaling thermostat to keep a constant 
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FIG. 9: The double-log plot of Cpp{r) for the volume fraction (j) = 0.093 with L = 112a and e = 1.0 
at i = 40ro as a function of the distance r. 
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FIG. 10: The double-log plot of Cpp{r) for the volume fraction (p = 0.37 with L = 72a and e = 1.0 
at t = 30to as a function of the distance r. 

temperature. Figures [9] and [10] are the results of the momentum correlation function of 
sheared elastic liquids with (p = 0.093 and (p = 0.37, respectively. These results also support 
the power law Cpp{r) ~ ^-s/s gyg^ ^^j^g elastic case, although Cpp{r) has the negative 
value for = 0.37 due to the back scattering effect. Although the existence of a power law 



Cpp{r) 



5/3 ]-^g_g i^QQYi believed 



28l | , it is the first time to verify the existence of such a tail 



from the molecular dynamics simulation. 
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VI. DISCUSSION AND CONCLUSION 



Now, let us discuss our result. We should discuss the validity of the generalized fluctuating 
hydrodynamics, (i) Can the noise be Gaussian? The answer is no. However, the non- 
Gaussian nature may not be important in our framework, because we only discuss the two- 
point correlation functions. When we discuss higher correlation functions, our theoretical 
treatment should be insufficient, (ii) What are contributions of the energy fluctuations? It 
might be important to include effects of the temperature fluctuation for denser cases. When 
we take into account such fluctuations, we must solve a cubic equation for the eigenvalues 
satisfying eq. (ITS!) to obtain the explicit form of Caik). The analytic expression of the 
solution is, however, rather complicated. Hence, we do not take into account the fluctuations 
of the kinetic temperature. This treatment may be justified as long as the USF is stable. 
Indeed, the temperature is immediately relaxed to be uniform, (iii) The contribution of 
inelastic collisions also appears through the inelastic Enskog operator for granular gases. At 
present, we have not solved the eigenvalue problem of the inelastic Enskog operator. In this 
sense, our treatment in this paper is far from the complete. The complete treatment will be 
discussed elsewhere. 

Serious contributions of inelastic collisions appear in large systems, because USF is un- 
stable for larger and highly inelastic systems. It should be noted that the shear flow is 
induced by the Lees-Edwards boundary condition. Using this boundary condition, USF is 
actually realized as shown in Fig. [TTl which is contrast to the shear flow obtained by mov- 
ing a physical wall [13]. However, our long-time simulation for larger systems is deviated 
from the theoretical prediction. In such a system, the linearized generalized hydrodynamics 
cannot be used, and the power law correlation function may disappear. This non- stationary 
tendency expected from unstable USF in larger systems can be observed in our simulation. 
Figure [12] shows the time evolution of Cpp{r) for the system sizes L = 36a and L = 72a in 
the case of e = 0.90 and (p = 0.37. It is clear that the result of a smaller system (L = 36a) 
converges, but the result of a larger system {L = 72a) does not converge. Although we have 



ts in 



51 



52 



arger systems. 



53|. 



not identified the reason why we could not obtain the converged resu 
it is likely to take place of the evolution of non-uniform structure [ ' 

It should be noted that the instability of USF can be reduced under the gravity, which is 
not included in this paper [sl, Q . Hence, we expect that the predicted long-range correlation 
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FIG. 11: The profile of the velocity u{y) along the flow direction as a function of the position y 
for cj) = 0.37, e = 0.90 and L = 32a with t = 20to. 



can be observed in experiments for sheared granular materials [18|, where the stable shear 
flow is realized. 

In this paper, we demonstrate that the framework based on the generalized hydrodynam- 
ics is useful to determine the structure factor or the pair-correlation function. This result 
is significant in the following two senses. First, we clarify the significant contribution of 
the shear to the structure factor of dense isothermal liquids, which is usually ignored in the 
papers to discuss glassy behaviors of sheared liquids [5J, l55|. Second, the determination of 
the structure factor gives complementary information to MCT for sheared granular liquids 

Q. 

In conclusion, we apply the generalized fluctuating hydrodynamics to isothermal sheared 
moderate dense liquids for both elastic cases and inelastic cases. The theory predicts that 
the density correlation function and the momentum correlation function, respectively, obey 
power laws r~^^^^ and r~^^^. Our theory can be valid even for the short-range scale for a < 
r < 5a. The density correlation function for the short range scale can be approximated by 
the theory introduced by Lutsko [6]. These results have been verified through the comparison 
with the simulation. 
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FIG. 12: (Color on line) The time dependence of Cpp{r) for the volume fraction 
t = 20ro , 40ro , 60ro for L = 36(7 and 12a in the case of e = 0.90. 
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APPENDIX A: APPROXIMATE EXPRESSION OF 5o(r,e) 

In this Appendix, we show the outline of unsheared pair-correlation function obtained 
by Lutsko Gj, who developed the generalized mean spherical approximation in describing a 
homogeneous cooling process of granular liquids. 

The approximate expression of (?o(e,T) is obtained from the inverse Laplace transform of 
the function G{t) satisfying 

Git) = / rfre-*V(7o(r,e), (Al) 
Jo 

where G{t) is given by 
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with the volume fraction = noO"^7r/6. Here, we introduce F{t) as 

F(f) = -{l + A,t + A,t') 

^' 12^{So + Sit + S2t^ + S3t^ + S^t^)' ^ ' 

where = 1, 5i = Ai - 1, S2 = - + 1/2, S3 = -A^ + Ai/2 - (1 + 20)/(120), 
54 = A2/2 - (1 + 2(p)/{l2(p)Ai + (2 + 0)/(240) with 



1 / (0-l)2-(60^o(a,e) + l)Z 
' 2^y 12<p{{2 + <p)-2go{a,e){<p-iyy 



A2 = goicT,e] 



(1 + 2(f))Ai - (2 + 0)/2 



(A4) 



1 + 6(j)go{a, e) 

In eq. (IA4p . the pair-correlation function (?o(o", e) at contact is given by the equilibrium pair- 
correlation function at contact Qeqio-) = (1 — 0/2)/(l — 0)^ as (7o(cr, e) = (1 + e)gcq{cr)/{2e). 
In addition, Z is given as 

-1 

(A5) 



where 



-1 



(A6) 



with P = ^(1 + 4(j)geq{(7))/m. 



APPENDIX B: TRANSFORMATION BETWEEN CARTESIAN COORDINATE 
AND OBLIQUE COORDINATE 

In the calculation of eqs. fl24|) - fH0l) . P7|) . and fl53l) . it is convenient to use the oblique 
coordinate, where the vector z of eq. in the Cartesian coordinate is transformed to the 
vector z in the oblique coordinate as 



where the matrix T is given by 



T 



z = Tz, 



1 

kx/k kykx/ {kk±) 
ky/k kj_/k 
kz/k —kyk^l {kk^) 



(Bl) 



kz/k± 


—kxj k\ 



(B2) 
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Here, we note that the second, the third and the fourth components of z are given by 
zi — ^3 = and ^4 = ^Usi where = ey • 5u, ^Ut — et ■ 5u, Sug — Ss ■ 5u, and 

u ^ 6u\\e\i + Sutet + SusCg (B3) 



with — {kx/k,ky/k,kz/k), ef — {—kykx/{kk±), k±/k, —kykz/{kk±)), and ej 
{kz/k±,0, -kx/k±). 

Prom this transformation, the right eigenvectors and the left eigenvectors in the obhque 
coordinate are respectively given by 

^^'^^ = ^(^fcr^o<J^A+,0,0), (B4) 

V'^'^^ = ^(^A;r^oa^A_,0,0), (B5) 

= (0,0,l,M(fc)), (B6) 

V'(')^ = (0,0,0,1), (B7) 

(B8) 

and 

= ^^,kp*{k,e),X+,Q,Q), (B9) 

= _L(,/,p*(/,,e),A_,0,0), (BIO) 
= (0,0,1,0), (Bll) 
= (0,0,-M(fc),l). (B12) 

(B13) 
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APPENDIX C: EXPLICIT REPRESENTATION OF C'„„(fc) 



In this Appendix, we explicitly derive eq. (IS^ . Cnnik) is related to Cnik) as Cnnik) 
a^^Cii{k). Using eq. fH5]) . we obtain 



Jo , _ 







F A+(fc(7*s))A_(fc(7*s))z/i*(A;(7*s))A;(7*s 



N+{k)N_{k) Ar+(fc(7*s))A^_(A;(7*s)) 
fc^ A+(A;(7*g))Vr(A:(rg))fe(7*g)^ 



N_{ky N_{k{rs))^ ■ ^ ^ 



From the relations 

^^-J^dsu*{k{r-s))k{rs)^ = -z/*(fc(7*t))^(7*t)2e~^o'^'"'^i(^(^*""^)^(^*'"^', (C2) 

^g-/o*rf^-A;{fc(r5)) = -A;(fc(7*t))e-^orf5A;(fc(7*5))^ ^(^3^ 

^g-/JdiAi{fc(ri)) ^ -X*_{k{'y*t})e-^o'i-s^UHrs))^ (C4) 

eq. fICip can be rewritten as 

r (J.-, ,,2,6f 4fc^A4fc)A„(fc) , k'X^{k)ul{k,e) , fc^A_(fc)z/*(fc,e) 



N4k)m4k)^ N+{k)^ N_{k) 

+An,a dt e 

X — 

dt 



4 



2 
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where we have introduced 

G'ii(^) - — j^^, — , m 

G22{k) - j;^, . (C8) 

It is possible to obtain the simpler expression for Cii{k) from the rela- 
tion {4eX+X_/{NlN^) + \+v{{k,e)k^/Nl + X^v{{k,e)kyNt) = {(T^nop*{k,e))-^ = 
So{k,e)/{o-^noA), and 



dt 

Ye-^^'''^''^^''''^'^^*-'^^^^^^^^^ (C9) 
dt { 

^*e-^«''^"^*T^«'''('(^*^"%n(M7*t)), (CIO) 
^*e-/o'^-"^T^^S^«*''('(^*^"))iy22(M7*^)), (Cll) 



with 



HMk) = k^^^ - + ^^'\k))GMk), (C12) 

Hn{k) = k^^^ - 2e\k)G^,ik), (C13) 
H22{k) = k^^f^ - 2e'\k)G,,{k). (C14) 
Substituting flCOD- dCTil) into eq. ([05]), we obtain 

C'n(fc) = noa^ 150(^,6) +7*Ai(fc)}, (C15) 

where we have introduced 

^ ' Jo N+ik)N4k) k{W 

Jo N+iky k{Yt)^ 

We, thus, obtain eq. (|52|) from eq. (IClSp with Cnn{k,t) = a^^Cii{k,i). 
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APPENDIX D: EXPLICIT REPRESENTATION OF Cpp{k) 

In this Appendix, we also explicitly derive eq. (1^ . The momentum correlation Cpp{k) 
can be written as 

Cpp{k) = {mn^a/tE?{C22{k) + C,,{k) + C,^{k)}, (Dl) 

where 

C22{k) + C,,{k) + C,,{k) = / diJ2 Yl i'^\k,t)^P^:-\-k,i)F^'"^\kiri)) 

"^0 l,m a=2,3,4 

E {-ir^'^^\k,t)^t\-k,t) 

l,m a=2,3,4 

xF('"^)(fc(7*f)), (D2) 
Here, to obtain ( 1D2I) we have used the relation 'if)^\k,i) =T-'0^'^(fc, t) and 

E„=2,3,4T«/.(fc)T„,(-fc) = (-l)^E„=2,3,4T^^(fc)T«.(fc) = 

Introducing 

POO 

D^{k) ^ / diJ2i^^\k,t)^^aH-k,t)F^^-\k{rt)), (D3) 

l,m 

we can rewrite Cpp{k) as 

C'pp(fc) = (mnoa/t£;)2{-D2(fc) + Ds{k) - D^{k)}. (D4) 
From ( ]D3[) . D2{k) can be rewritten as 



X 



iV+(A;)A^^(A;) A^+(A;(7*s))A^_(A;(7*s)) 
+2Ay""rftV/o'^^-n(Mrs-))e-^''^^"^*Ti^«'''('^(^*^")) 



X 



+2A I die-fod''-^Hr-s))e-I^'^^r'-^^i^ 







A-(A:)^ A-(fc(r^-))Vi(fc(r^-))A:(r^-)^ 
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From the parallel procedure to that in Appendix [Cl we obtain 

- ^\N^{kfN^{ky^ N^{kY^ N.{kY] (D6) 

where 

^ ~ Jo N^{k)N_ik) kirt)- 

Jo N^iky km' 

r Xi{krH22Ckm) k^~ky{rt) , fu^x^^ukr^)) 
^Jo N_{kf km' ■ ^ ^ 

From AX+{kf\^{kf/{N+{kfN_{kf) + \+{kY/N+{kf + \_{kY/N_{kf = 1, we find that 
D2{k) is represented by eq. flM|) . 
From (]D3p . we obtain 

D,Xk) = 2A dte-^^'^~''''*^^''^^*'^^^^^*'^e^°^^^^^^v*{k{^^ (D8) 
Jo 

From the identities 



and 



e/^-*^ = (DIO) 



we find that D3(k) is expressed as eq. fl55|) with 



oo 



k^ 



From eq. fIDSp . we obtain 

L)4(fc) = -2A |M^(fc) ^^^y - 2M(fc)M(fc(7*g))^^|^ + M(fc(7*g))^ + 1 

^^-J^ds2uUHrmrs)\*^^Yt))km'. (D12) 
With the aid of flD9p . we confirm that Di[k) can be represented by eq. fl56l) . where 

A4(/c) = rft I ^-.^^(^^^'^ F(fc, t) - \ F{k, t)e- ^o^'""^ (Hrmir^? (D13) 

^ ^ Jo \ km' km J 

and t) = M{k{t)) - k{t)M{k) jk. 

Finally, substituting eqs. (15^ . (!55|) and (156|) into (IDip with (lD2p . we obtain the explicit 

representation of Cpp{k) as eq. (!53l) . 



27 



APPENDIX E: EVALUATION OF ASYMPTOTIC FORM OF Aj{r) 



In this Appendix, we evaluate the asymptotic form of Aj(r) for r ^ Ic- 

1. The Asymptotic Form of Ai(r) 

We note that Ai(r) can be separated into three parts 

Ai(r) = An(r) + Ai2(r)+Ai3(r), (El) 

where 

^ ^ r dk sm{kr/a) Ak''H^2{km))kMi*t) 
"^"^^ n,a^J (27r)3 kr/a J, N+{k)N_{k) k{rt)' 

X e" /o (fc(7*5))+A(2) (fc(r 5))} ^ (E2) 

A - f sin (kr/a) ^- PHniHri)) k^'kyj^t) _^ dsxW ,)) ^^.o^ 

- n,a^ J (27r)3 kr/a J, N^{ky k{YiY ' 

A M - ^ / dk sin(fcr/a) r ^ k^H22{Hrt)) k,~ky{ri) ^^2 fUsX^'UMrs)) /p.N 

Let us introduce Ic = (7/^/7*, = r/lc, and the transformation i = r/7*, s = t'/7*, 
fc = \p)*K.jr*' = aK/r, and fc(r) = y/^K{T)/r*. By using the new variables, eqs. (]E2p . 
(IE3p . and (IE4p can be rewritten as 

^^^^^^ ^ ^*-5_^ f dK_s^ r AKm,2iV¥KiT)/r*) K^Kyir) 



2 



^g-2/o-dr'r-^A(i)(Vr^(-')A*)^ (E6) 



^ r*-^_^ r^sin^ r^^K'H22{V¥K{T)/r*)K^Ky{T 



noaHl J (27r)3 Jo A^_(v^fs:/r*)2 i^(r)2 

Introducing the transformations t = r 

*2/3y and r' = 7-*2/3y/ extract the asymptotic 



28 



forms for r* ^ 1, Aij(r) (j = 1, 2, 3) can be rewritten as 
A 



Aii(r) = r 



*-13/3_ 



dK sinK r AK^Hi2{V¥K{r*^/^Y)/r*) K^Ky{r*^/'^Y) 



^g-/;dy'r-V*2/3{A(i)(v/^K(r*2/3y')/^.)+A(2)(Vr^(r*2/3y')/r*)}^ (E8) 



Ai2fr) = r 



*-13/3_ 



A 



A 



(E9) 



r 



*-13/3_ 



dK sinK ^^K^H22{^K{r*^I^Y) /r*) K^Ky{r*^I^Y) 



noaHl J (27r)3 K Jq ^ N_{^K/r*Y 

^ g-2 /J" dy'7*-lr*2/3 A(2) {v^fc(r*2/3 Y')/r* ) 



(ElO) 



where the functions included in (lESp , (1E9P , and (lElOp have the asymptotic forms for r* ^ 1 



Ky{f^/^Y) - r*2/3{y7^^ + 0(r*-2/3)}, 



i7ii(v/l^ir(r*2/^r)/r*) 
i/22(Vri^(r*'/='r)/r*) 



2co 
2^5^ 



y\K,\r*-'/' + 0{r*''/'), 



(Ell) 
(E12) 

(E13) 
(E14) 
(E15) 

(E16) 
(E17) 



■ /J- dy'7*-ir-2/3A(l)(VrK(r-2/3y')/r.) _ ^^07*-^/^ /o'' dy'\K,\Y'r*^/^ 



K 



xe 



|ir^|r*2/3y 

-/„^dy'.j(o,t)K^y'2 1^ ^ o(r*-2/3)} , (E18) 



- dy'7'-lr*2/3A{2){Vr ^ pico7*-^/2 dY'\K^\Y'r*^/^ 



K 



1^ |^*2/3y 



^^-!^dY'uimKlY'2 1^ ^ o(r*-2/3)} , (E19) 
with Cg = nocr^p*(0, e). Substituting these relations into eqs. flESp - flElOp . we obtain 

An(r) = Jnr*-i^/3^0(r*-i3/3), 
|Ai2(r)| = |Ai3(r)| = Jisr*-^ + 0(r*-i^/3), (E20) 
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where we have introduced 



2; I 

00 



4noa3/3 y^c^y (27r)3 F 
Thus, the asymptotic form of Ai(r) can be represented by eq. fl60|) . 

2. The Asymptotic Form of A2(r) 

Similarly, A2(r) is given by summation of the three terms 

A2(r) = A2i(r) + A22(r) + A23(r), (E23) 

where 

A2i(r) = / sin(fcr/a) j-^ ^-AX4k)X4k)HuCHrt}) KkyiTt) 



(27r)3 A:r/(T iV+(A:)Ar_(A;) ^(7*t)2 

xe-i'o'^""^^*''(^(^*^"))+^'"(^(^*""))>, (E24) 
. _^ I dk sm{kr/a) ,-X+{kf HiiCHYt}) k,~ky{ri) 

^221^") — O" / — — : — ; / M 



(27r)3 kr/a Jo N+{k)^ ~k{Yt)^ 

^^-2j^dsxWCHr-s))^ (E25) 

^ f dk sm{kr/a) r ,-Xi{kf H22CHrt}) k,~ky{rt} 
^23 — O" / TZ — TT- : ; / dt 

xe" 



(27r)3 kr/a Jo N.^kf k{'y*t}^ 



As in the case of the previous subsection, from the transformations t = r*'^/^Y / 'j* , s 
r*^/^Y'/Y, and k = ^/7*K/r*, we obtain 

^ _ r*-^/3 gjj^^ ^^AX+{y^K/r*)X4y^K/r*)Hi2iV¥Kir*^/^Y)/r*) 



Yll J (27r)3 K Jo N+{^K/r*)N_{y/YK/r*) 



dY'r-^r*y'^{XW(^Kir*2/'iY')/r*)+X(2\^k{r*y'-'Y')/r*)} (E27) 



r*-^/3 r gi^^ /-"^ A+(v^ir/r*)2i7n(v^i^(r*2/3y)/r*)ir^iry((r*2/3y) 

A22('") — ~ TTT- I — — / dV- 



7*/3 J (27r)3 K Jo N+{y/YK/r*y ^((^*2/3y)2 

^ g-2 dy 7- V*2/3a(i) (VrK(r*2/3yO/r*) ^ (^£28) 

- r*-^/3 /. gjj^^ p A_(v^K/r*)2if22(v^i^(r*2/3r)/r*) ir^.iry(r*2/3y) 



7*/3 _/ (27r)3 K Jo N_{y/YK/r*y K{r*^/^Y) 

^g-2/J'dy'7*-iA(2)r*2/3(^K(r*2/3y')/r*)_ (E29) 
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From eqs. flETTl) - flET9|l and eqs. flE27D - flE29|) . we obtain 



A2i(r) = ^^Ji,r*-"/3 + 0(r*-i3/3), (E30) 
|A22(r)| = |A23(r)| = "-^J^.r*-' + 0{r*-''/% (E31) 
We, thus, find the asymptotic form of A2(r) is given by eq. (!6T!) . 

3. Asymptotic Form of A.^{r) 
A3(r) can be represented by 

J {211^ kr/a Jq 



From the transformations t = r/j*, s = r'/'j*, and k = y^K/r, we can rewritten ( 1E32I1 



as 



Introducing the other transformations r = r*'^/^y and r' = r*'^^^y\ we obtain 



(27r)3 K ^ 



xe ^ . (E34) 

We should note that 



Z/n ■ ~ U, 



(0)+O(r*-^/=^) (E35) 



can be used for r* ^ 1 from eq. flElip . Using this equation and eq. flE12p . we obtain 



As{r) = J3r*-5/=^ + 0(r*-2/3), (E36) 

where 

Therefore, the asymptotic form of A3(r) is given by eq. fl5^ . 
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4. The Asymptotic Form of A4(r) 



AJr) can be rewritten as 



A4(fc) = a-' [ . "^^^3 (^^/^) r dm{k,i)e'^o'i''''*-iHrmir-s)\ (ess) 

J (ztt) kr/a jq 

where we have introduced F/^{k^t) = |^||^^^^F(A;, t) — -^^^^^ | F(k,i). From the transfor- 
mations i = r*'^/^Y/Y, s = r*'^/^Y'/Y, and k = ^/YK/r*, we obtain 

Substituting the asymptotic expression 

F4(VrK/r*,r*2/3y/^*) ~ I!!^ + i^M(K)^ r*^/^Y (E40) 

\ 2 Ax A J 

for r* > 1, eqs. flETTD and flET2D into eq. flE39D . we obtain 

A4(r) = J4r*-^/=^ + 0(r*-2/3)^ (^4^^ 

where 

Hence, we find that the asymptotic form of A^lr) is given by eq. fl63p . 



[1] J-P. Hansen and McDonald, Theory of Simple Liquids 2nd Edition (Academic Press, London, 
1986). 

[2] U. Balucani and M. Zoppi, Dynamics of the Liquid State (Oxford University Press, Ox- 
ford, 1994). 

[3] I. S. Aronson and Lev S. Tsimring, Rev. Mod. Phys. 78, 641 (2006). 
[4] O. Pouliquen, Phys. Fluids 11, 542 (1999). 

[5] L.E. Silbert, D. Erta§, G. S. Grest, T. C. Halsey, D. Levine and S. J. Plimpton, Phys. Rev. E 

64, 051302 (2001). 
[6] J. F. Lutsko, Phys. Rev. E 63, 061211 (2001). 

32 



[7] M. Alam and S. Luding, Phys. Fluids 15, 2298 (2003). 

[8] GDRMiDi, Eur. Phys. J. E 14, 341 (2004). 

[9] J. F. Lutsko, Phys. Rev. E 70, 061101 (2004). 

[10] N. Mitarai and H. Nakanishi, Phys. Rev. Lett. 94, 128001 (2005). 

[11] V. Kumaran, Phys. Rev. Lett. 96, 258002 (2006). 

[12] A. V. Orpe and A. KudrolH, Phys. Rev. Lett. 98, 238001 (2007). 

[13] K. Saitoh and H. Hayakawa, Phys. Rev. E 75, 021302 (2007). 

[14] N. Mitarai and H. Nakanishi, Phys. Rev. E 75, 031305 (2007). 

[15] H. Hayakawa and M. Otsuki, Phys. Rev. E 76, 051304 (2007) . 

[16] M. Otsuki and H. Hayakawa, larXiv:0711.T42T] 

[17] H. Hayakawa and M. Otsuki, Prog. Theor. Phys. 119, 381 (2008). 

[18] A. V. Orpe, V. Kumaran, K. A. Reddy, and A. Kudrohi. [aFXiv:0811.2848] (2008). 

[19] H. Jaeger, S. R. Nagel, and R. Behringer, Rev. Mod. Phys. 68 , 1296 (1996). 

[20] N. Brilliantov and T. Poschel, Kinetic Theory of Granular Gases (Oxford University Press, 
Oxford, 2004). 

[21] J. T. Jenkins and M. W. Richman, Phys. Fluids 28 3485 (1985). 

[22] J. T. Jenkins and M. W. Richman, Arch. Ration. Mech. Anal. 87, 355 (1985). 

[23] V. Garzo and J. W. Dufty, Phys. Rev. E 59, 5895 (1999). 

[24] J. F. Lutsko, Phys. Rev. E 72, 021306 (2005). 

[25] H. Hayakawa and M. Otsuki. larXiv:0805.1515l (to be published in Prog. Theor. Phys. Suppl.) 

[26] M. Otsuki and H. Havakawa. larXiv:0805.1516l (to be published in Prog. Theor. Phys. Suppl.) 

[27] V. Kumaran, to be published in Phys. Rev. E (Dynamics of a dilute sheared inelastic fluid. I 
& n) (2008). 

[28] J. F. Lutsko and J. W. Dufty, Phys. Rev. A, 32, 3040 (1985). 

[29] J. F. Lutsko and J. W. Dufty, Phys. Rev. E, 66, 041206 (2002). 

[30] H. Wada and S. L Sasa, Phys. Rev. E, 67, 065302(R) (2003). 

[31] B. M. Law, R. W. Gammon, and J. V. Sengers, Phys. Rev. Lett. 60, 1554 (1988). 

[32] P. N. Segre, R. W. Gammon, J. V. Sengers, and B. M. Law, Phys. Rev. A 45, 714 (1992). 

[33] P. N. Segre, R. W. Gammon, and J. V. Sengers, Phys. Rev. E 47, 1026 (1993). 

[34] J. R. Dorfman, T. R. Kirkpatrick and V. Sengers, Ann. Rev. Fluid Mech. 45, 213 (1994). 

[35] W. B. Li, J. V. Sengers. R. W. Gammon, and P. N. Segre, Int. J. Thermophys. 16, 23 (1995). 

33 



[36] W. B. Li, K. J. Zhang, J. V. Sengers, R. W. Gammon, and J. M. Ortiz de Zarate, Phys. Rev. 
Lett. 81, 5580 (1998). 

[37] W. B. Li, K. J. Zhang, J. V. Sengers. R. W. Gammon, and J. M. Ortiz de Zarate, J. Chem. 
Phys. 112, 9139 (2000). 

[38] M. Otsuki and H. Hayakawa, to be published in Proceedings of Rarefied Gas Dynamics 2008 

(|arXiv:0g09.1118p . 
[39] J. F. Lutsko, Phys. Rev. E 66 , 051109 (2002). 

[40] W. W. Gotze, in Liquids, Freezing and Glass Transition, edited by J. P. Hansen, D. Levesque, 

and J. Zinn-Justin (North- Holland, Amsterdam, 1991). 
[41] S. P. Das, Rev. Mod. Phys. 76, 785 (2004). 
[42] D. R. Reichman, J. Stat. Mech. P05013 (2005). 
[43] T. R. Kirkpatrick, Phys. Rev. A, 32, 3130 (1985) 

[44] T. R. Kirkpatrick and J. C. Nieuwoudt, Phys. Rev. A 33, 2651 (1986). 

[45] S. P. Das, Phys. Rev. A 42, 6116 (1990). 

[46] S. Sinha and M. C. Marchetti, Phys. Rev. A 46, 4942 (1992). 

[47] I. M. de Schepper and E. G. D. Cohen, J. Stat. Phys. 27, 223 (1982). 

[48] W. E. Alley and B. J. Alder, Phys. Rev. A 27. 3158 (1983). 

[49] T. P C. van Noije, M. H. Ernst, R. Brito, and J. A. G. Orza, Phys. Rev. Lett., 79, 411-414 
(1997) 

[50] M. Alam, P. Shukla, and S. Luding, J. Fluid Mech. 615, 293, (2008). 

[51] M. Alam and P. R. Nott, J. Fluid Mech. 377, 99, (1998). 

[52] M. Alam and P. R. Nott, J. Fluid Mech. 343, 267, (1997). 

[53] M. Lee and J. W. Dufty, Phys. Rev. E 56, 1733 (1997). 

[54] M. Fuchs and M. E. Gates, Phys. Rev. Lett. 89, 248304 (2002). 

[55] K. Miyazaki, D. R. Reichman and R. Yamamoto, Phys. Rev. E 70, 011501 (2004). 



34 



